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Abstract. Along with data on radial velocities, more and more 
data on proper motions of individual stars of globular clusters 
are becoming available. 

Their usage was until now rather limited. It was mostly re¬ 
stricted to determining cluster membership of individual stars 
and to determining the spatial velocity of the cluster. We will 
study the two dimensional distribution of the proper motions, 
in order to clarify the relation between the projection of the ve¬ 
locities and the dynamical structure. Obviously some assump¬ 
tions on the dynamics of the system have to be made. In this 
study we chose a Plummer model. This complements an earlier 
study about line profiles. 

The three velocity dispersions in the main directions are de¬ 
termined and compared with the three spherical velocity dis¬ 
persions. Two parameters analogous to Binney’s anisotropic 
/i-parameter are defined. Moreover, since the proper motion 
components are distance dependent and the line of sight ve¬ 
locity dispersion is distance independent, the distance can be 
determined. 

Key words: Celestial mechanics, stellar dynamics - globular 
clusters: general 


1. Introduction 

One of the greatest challenges in astrometry is the accurate 
measurement of minute relative variations in the positions of 
objects in the sky, in order to measure proper motions. Recent 
advances in data acquisition techniques, largely symbolized by 
the Hipparcos Mission, will produce a large quantity of proper 
motion data of high quality. The more traditional approach 
simply takes advantage of the lapse of time: the longer the 
time span that is covered, the better proper motions are de¬ 
termined, at least if the first-epoch medium is sufficiently well 
behaved. No wonder that useful proper motions are also be¬ 
coming available in those areas of the sky that happened to be 
of particular interest to astronomers (then and now). Globular 
clusters are such areas, and plates exist that are almost a cen¬ 
tury old. There is already a body of data on proper motions 
for globular clusters (Cudworth fc Smetanka 1992, Rees 1992 
and references therein). Similarly, other data sets, located in 
popular fields in our own Galaxy and the Galactic Bulge, are 
being produced (Minniti 1992). 

In the case of globulars, one of the uses of proper motions 
consists in the determination of cluster membership by elimi¬ 


nating stars with anomalously large proper motions. From the 
astrophysical point of view, such a procedure produces cleaner 
color-magnitude diagrams. Another more recent use is the de¬ 
termination of the proper motion of the globular itself (Cud- 
worth and Hanson 1993). On the other hand, proper motions 
are useful for what they directly mean: they provide kinemat- 
ical information. They are similar to the line-of-sight veloc¬ 
ity, which has proven to be very useful in galactic cluster and 
galaxy research. Therefore we set out to study the theoretical 
aspects involved in calculating and interpreting proper mo¬ 
tions. This we do very much in the same way as was done by 
one of us for the line-of-sight velocities (Dejonghe 1987, here¬ 
after Paper I). 

In sections 2 and 3, we consider a one-parameter sequence 
of Plummer models, which covers a variety in orbital struc¬ 
ture. It is well-known that the Plummer model is not fit to 
describe stellar systems, be they globulars or spherical galax¬ 
ies. Nevertheless, we consider the Plummer model because of 
its analytical simplicity, and we are confident that our results 
are applicable to real stellar systems, at least in a qualitative 
sense. A similar statement can, in hindsight, be made about 
the models in Paper I. 

We study the projected proper motions and compare our 
results with the earlier calculated line profiles for each member 
of the Plummer family. Conversely, we show what dynamical 
information can be extracted from observed proper motions. 
In particular, it was our goal to get information on the or¬ 
bital structure (i.e. radial or tangential orbits) of the system, 
regardless of our distance to the system. 

Since, as was already mentioned, the distance-dependent 
proper motions are very similar in nature to the distance- 
independent line-of-sight velocities, it is clear that the simul¬ 
taneous knowledge of line-of-sight velocity distributions and 
proper motion distributions will constrain the distance to the 
system (section 4). A first attempt at this can be found in 
Lupton, Gunn fc Griffin (1987). 

Much of the work can be carried out analytically, which is 
not only a useful check on the numerical calculations, but gives 
this work, we hope, an additional quality which is increasingly 
rarely found. 

2. The theoretical framework 

2.1. The Plummer family 

We consider anisotropic distribution functions F(E,L ), which 
are functions of the specific binding energy E = 0( r) — ^( v^. + 
v ii + rj,), whereby i/'(r) = —V(r) > 0 is the potential function 
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the angular momentum vector. 

The Plummer model has the well known potential 


4'(r) = (1 + r 2 


( 1 ) 


expressed in some suitable units, with the corresponding mass 
density 


P(r) = ^(l + r 2 ) l (2) 

Dynamical anisotropic models for the Plummer mass den¬ 
sity can be constructed considering the augmented mass den¬ 
sities 


p{f,r) = -^-f b q {l + r 2 ) = , 


(3) 


dependent on the parameter q , but yielding the same physical 
mass density for all q , which is obtained by substituting Eq.(l) 
in Eq.(3). The distribution function F(E,L) that corresponds 
to a p(4\ r) turns out to be (Dejonghe, 1986) 
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L 2 > 2 E (5) 


where Tfi) denotes the gamma function. Positiveness of 
F(E,L) restricts the value of q to q < 2. The physical mean¬ 
ing of q is most obvious when considering Binney’s anisotropy 
parameter (Binney fc Tremaine, 1987, p204) 
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The important point to note is that q and j3 have the same sign. 
Hence radial orbits are prevalent if j3 > 0, and more tangential 
orbits are prevalent in a /3 < 0 cluster. If-jf = 0 then both types 
occur in equal amounts in an isotropic sense. If q becomes very 
negative, then cr~ — r 0, so in the limit q —*■ —oo the stellar 
system only contains circular orbits. The other limiting case 
q = 2 has a /3 which tends to 1 for great distances of r, so in 
the outer regions only radial orbits are present. We classify the 
models into radial systems, with 0 < q < 2, tangential systems, 
with q < 0, and isotropic systems, with q = /3 = 0. The latter 
case is normally called a Plummer model. 


2.2. Notations - normalisations 

The velocity v of a star in spherical coordinates is in obvious 
notations: v = v r e r + Vpep + v q )et). Generally we cannot observe 
these components: at best we obtain the projection of v along 
the line of sight which we denote by v p and the proper motion 
which is the projected velocity on the plane of the sky. This 
proper motion is normally measured in arcsec/year in some 
coordinate system, say (p a . Fs)- For this paper it is useful to 
define a system of perpendicular axes so that the proper motion 
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projected radius, and pt which is measured perpendicular to 
it. The velocity components v c and Vt are obtained by the 
relations v c = p c A and vt = Pt A, whereby A is the distance 
to the system. In the following, we will refer to v c as the center 
velocity (and not radial velocity, since otherwise we would add 
to the confusion about the term radial) and vt as the transverse 
velocity. So in this system the velocity v is: v = e c e c + vtGt + 
v p e p . We will transform the frame (e r , e v , ) to (e c ,et,e p ), 
considering the isotropy in & and ip. The axes e r ,e c and e p 
form a plane a. We rotate the axes and e P so that e P 
lies in a. Then et is lying according to e t ?, since both vectors 
are perpendicular to a. Let r be the distance and r p be the 
projected distance from the center of the cluster to the star; 
and let s be the algebraic distance from the star to the plane 
of the sky which cuts through the center of the system. Figure 
1 gives a perspective view of both coordinate frames and the 
distances r, r p and s. The transformation is a simple rotation 
around et = e,j and the formulae are then: 
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Fig. 1. The geometrical relations between the frame related to the 
spherical coordinate system (e^e^e,?) and the frame related to 
the plane of sky and the line of sight (e c , Gt,G p ). 

The projection f p (r p ) of any scalar quantity /(r) on the 
plane of the sky is 

/ +oo /*+CO J 

f(r)dz = 2 f(r) If— . (8) 
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For example, the projected mass density in the Plummer model 
reads: 

Pp(t p ) = ^(1 + r 2 v )~ 2 = (9) 

7T 7T 

whereby 4’p = 4’{' r p) = (1 + r p )~^ is the potential on the 
projected distance r p from the center. 

Hence, the cumulative mass as a function of r p is 


M p (r p ) = 2n J Pp{rp)r p dr p = 


1 + r 2 p 


( 10 ) 
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which can be used to define the core radius as the radius at 
which p p (r p ) drops to half of its central value: r 2 c = \f2 — 1 —»■ 
r c ~ 0.64 . 

To determine the proper motion distribution, it is conve¬ 
nient to redefine the independent variables (the projected ve¬ 
locities): at every radius, the velocity is normalized to unity 
for the largest proper motion corresponding to bound orbits. 
The maximum velocity which occurs at the distance r p from 
the center is given by u max = \/2 ip p . The normalized velocity 
components x, Vt C are thus 
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3. Proper motion distributions 


pm r (v c , v t ) dv c dv t 


= dv c dvt 


dr 2 


V r2 - 


F(E, E) dv p 


need not be specified by the double prescription (4) and (5). 
The PM-distribution (15) becomes after normalisation: 
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where / max is the value of / for which E equals zero, or / max = 
((tp/tpp) — v 2 ) 1 / 2 . The integration over the velocity component 
/ results in 
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and we also denote the normalized modulus of the proper mo¬ 
tion by 
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Note also that this normalization dispenses with the need to 
distinguish between the proper motion (p c , p-t) and the proper 
motion velocities (y c , vt). The energy and angular momentum 
expressed in the normalized velocities are: 


To avoid cumbersome notation, we set h = 2ra — m + j + 4. For 
the integration through the cluster we need to calculate 
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With the substitution x = ip p { 1 + r 2 ) 2 this leads to 


(14) 
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, l' (±-.’)V-*"-‘^. ( 20 ) 

The powers of the last two factors are positive integers so they 
can be expanded. Upon interchanging the integration and the 
summation signs, and using Euler’s integral the above equation 
transforms into: 

, 2) - 2m - 1 r(»-j + |) (i E'~ 3+i 


Let pm rp (» c ,»() be the proper motion distribution (PM- 
distribution) at r p . To calculate this PM-distribution we need 
to perform a double integration. The first integral averages over 
the velocity-component along the line of sight (v p ) and the sec¬ 
ond integrates through the cluster (z). The PM-distribution is 
given by: 


E = ip : 


(15) 
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and depends of course on the distribution function F(E,E). 
The boundary r max is the value of r for which the escape ve¬ 
locity is reached, or: 1 + 70 ax = />“ ax = 1 /( 0 ; 0 ). 

ft is useful to normalize the proper motion probability dis¬ 
tribution to one, though it really is proportional to the pro¬ 
jected mass density given in Eq. (9). 
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The PM-distribution (17), using polar coordinates namely 
X = v cos a and 77 = v sin a and a transformation of the hy¬ 
pergeometric function then becomes: 


3.1. Proper motion in the case q = —2 n (n a natural number) 

In this case the PM-distribution can be calculated exactly. The 
hypergeometric series in the distribution function F(E, E) is 
then given by 
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This result can be made considerably simpler for zero velocity. 
In this case integration over the velocity component / as given 
in Eq. (18) yields: 
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since all positive powers of g and x vanish. The integration 
through the cluster results in a beta-function, so that the PM- 
distribution for zero velocity is given by: 


(P m r p (M))q=-2n 
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This expression can also be obtained from Eq.(22) by taking 
the limit for v —»■ 0. 


Equation (22) can also be expressed in elementary func¬ 
tions for specific q, using Gradshteyn and Ryzhik (1965, Eqs. 
9.121). The simplest case is the isotropic one (q = 0), which 
reads: 


Pm rp (x,»?) 
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If v —»■ 1, then pia r ^(x, v) goes to the expected value zero. If 
v —»■ 0, then pm r (0,0) = 15/8. Equation (25) shows that in 
this case the PM-distribution does only depend on the velocity 
v 2 , which is not surprising. More surprising however is that all 
normalized PM-distributions are the same irrespective of the 
projected distance r p . 

This result holds also for other distribution functions 
within the same potential, but which, of course, do not pro¬ 
duce a self-consistent mass density. For example if F(E) = 
A(2E)~ + m , with A some constant, the integration I c through 
the cluster and over the velocity component £ (we work again 
with normalised velocities) results in: 
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The spatial mass density is given by: 


p(r) = 47r I F(E) v 2 dv = 47ry4(2^) m " 1 " 2 i?(— , ra + —) 


(27) 
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P P (r p ) = 87T+(2/> p ) m+1 B(y m + ^)B(i, y + i). (28) 

One easily checks that for m = 3, using the appropriate analyt¬ 
ical expressions of the hypergeometric function (see Gradsteyn 
& Ryzhik, 1965), the original isotropic Plummer model is re¬ 
covered. 


Another particular case is q = —2. Obviously, this is a 
more elaborate example. The PM-distribution (22) expressed 
in elementary functions, leads after some calculations to: 
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Again, if v —»■ 1, then pm rp (y, i/) goes to the expected value 
zero. And if v —»■ 0, then pm r (0, 0) = 63(l + />/)/64. The result 
here depends on the projected distance and on the velocity 
components. Later we will use this exact distribution to test 
our numerical results. 

This kind of calculations could not be performed with the 
line profiles (see paper I), since for q = —2 n the power of 
Cmax in Eq. (18) is not an integer so this expression cannot be 
expanded. 


3.2. Proper motion in the general anisotropic case 


The distribution function F(E, L ) is given by Eqs. (4) and (5). 
The integral J F(E, L) d( cannot be calculated in a general 
manner due to the cumbersome expression of F(E, L ) as a 
function of £, so we have to work with approximations. How¬ 
ever, we can calculate exactly 


1. the “top value”: this is the probability of Ending a star at 
r p without proper motion, i.e. pm rp (0, 0), and this is found 
to be (see appendix A): 
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If q = —2n this equation reduces, after a transformation of 
the hypergeometric function, to Eq.(24). 

2. all the moments of the two-dimensional PM-distribution. 
They are determined in appendix B and are: 
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whereby the normalized moments are given by 
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Fig. 2. Probability contours for the PM-distribution of different 
radial g-clusters. Left panels at the core radius ( r p = 0.6), right 
panels at five times the core radius. 


Fig. 3. Probability contours for the PM-distribution of different 
tangential g-clusters. Left panels at the core radius ( r p = 0.6), right 
panels at five times the core radius. 


In particular, two of them have special significance in view 
of later applications: 
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It is a logical step to try to use the above results as an aid 
to calculate the PM-distributions. To that end, we conjecture 
that the distribution can be well approximated by the series 


P m r p (\ i *7) = ^2 aij (1 " ^ " , ? 2 ) a+ ' *? 2J - (35) 

’<3 


We can determine the a, 3 coefficients by matching the moments 
of the approximated distribution in Eq. (35) with the true mo¬ 
ments and by equating pm (0,0) from Eq. (30) to y~7 a, p. 
Therefore it is useful to obtain an analytical estimate for the 
exponent a. 

For small values of r p the PM-distribution is approximated 
by (see appendix C) : 
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and for large values of r p this approximation is: 
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In Figs. 2-3 we plot the contour lines of the distribution 
for some values of q and r p . The calculations of the two di¬ 
mensional PM-distribution were checked twice, once by nu¬ 
merical integration at some points and once by comparison of 
the results with the exactly calculated distributions for q = 0 
(isotropic case) and for q = — 2. The matches were excellent for 
tangential systems as long as the q values were not too negative 
{q < —6) and the distance r p was not too large. In the other 
cases, the individual terms in the calculation of the moments 
(32) tend to become very large, causing subtractive cancella¬ 
tion, so that the numerical solution for the coefficients and the 
subsequent expression (35) is correspondingly unreliable. For 
radial systems the series approximation did not deliver such 
good results for large positive (> 1) q values or large r p values. 
The reason is that the PM-distribution does not approach zero 
uniformly on circles. Hence the figures that did not give a good 
match due to accuracy errors were produced by straightforward 
(but fairly slow) numerical integration. 

Admittedly, the surfaces in Figs. 2-3 are only exact for the 
Plummer models we consider here. It is not unreasonable to 
assume, however, that our results carry over, at least in a qual¬ 
itative sense, to more realistic models. 

We immediately notice the difference between the PM- 
distributions for both cluster types. For the radial systems 
(Fig.2) the PM-distribution is unimodal and becomes more 
and more elongated along the center-axis. This feature is more 
prominent with increasing q and r p . This can be explained by 
the fact that our radial clusters have predominantly radial or¬ 
bits in their outer parts. There is no contamination of the (spa¬ 
tial) tangential velocities in the center direction. For tangential 
systems, (Fig.3) the PM-distribution becomes bimodal towards 
the transverse direction. This feature is again more prominent 
with increasing \q\ and r p . Again, this can be explained by not¬ 
ing that our tangential clusters have predominantly tangential 
orbits in their outer parts. There is no contamination of the 


Hence we adopt a = 9/2 — q. 














Fig. 4. These pictures show different cuts of the PM-distributions along both velocity axes, at the core radius (r p = 0.6, left panels) and 
at five times the core radius (right panels). Values of q are indicated in the figure. 


(spatial) radial velocities in the tangential direction. 

It is important to realise that Figs. 2 and 3 provide an im¬ 
mediate way of determining qualitatively the orbital structure 
of star clusters! There is no need for a distance determination. 
In particular we note the striking difference between the right 
panels of Figs. 2 and 3, which show distribution at about 5 
core radii, which is therefore a radius where a clearcut diagno¬ 
sis should be possible if sufficient stars are present. A different 
pair is supposedly the lower-left panel of Fig. 2 and the upper- 
left one of 3 showing the PM-distribution at core radius of two 
clusters with not too different an orbital structure. Neverthe¬ 
less, it does not seem hopeless to distinguish both with good 
quality data. Clearly the dispersions of the marginal distribu¬ 
tions must play a key role, and we discuss them in the next 
section. 

Finally, Fig. 4 shows, more quantitatively, the cuts of the 
distributions along both velocity axes. Here we see the depen¬ 
dence of the PM-distribution upon \ or r] for various q. The 


more positive q , the more peaked the PM-distribution is in 
the top. The more negative q , the more bimodal the distribu¬ 
tion becomes. The effect is more pronounced at large projected 
radii. 


4. The dispersions 

The second order moments in the proper motions are of partic¬ 
ular interest, since these are the most securely determined. The 
first order moments in these models are zero. Using Eq. (32), 
we find for the dispersions 
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Fig. 5. The line of sight velocity dispersions and the transverse velocity dispersions are plotted for various q. Note the qualitative similarity. 
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This last result is very similar to the line of sight velocity dis¬ 
persion calculated in paper I: 
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In Fig.5 the line of sight velocity dispersion and the trans¬ 
verse dispersion are plotted against the radial distance r p , for 
varying q. Both have qualitatively a very similar behaviour. 
In both panels of Fig.5 there are fixed points, due to the fact 
that q enters linearly in both the numerator and the denom¬ 
inator of Eq. (39) and Eq.(40). For the line-of-sight velocity 
dispersions the fixed point is reached at r p = y/3/2 and has 
the functional value a p = (3tt/ 64) y/3/5. In the transverse 
direction the fixed point is located at r p = 1/2, with func¬ 
tional value of = 3 7r/ (32 \/E). Also a maximum occurs for <r p 
at r p = —(12 + 10g)/(12 — 5g) if g < —6/5 and for <r/ at 
r p = —(12 + 9g)/(12 — 6g) if g < —4/3. Both figures have in 
common that radial clusters look hotter at the center than tan¬ 
gential clusters, in both projections. On the contrary, at the 
edge the tangential clusters look hotter, as can be expected. 
The difference lies in the magnitude of the effect, which is 
greater for cr p than for a"j: in the center, while the opposite is 
true at larger distances. 

The center velocity dispersion on the contrary is a simple 
monotonous decreasing function of r p . Moreover, for every r p 
it is an increasing function of the parameter q. So at each 
projected distance radial clusters give a higher contribution to 
the center velocity dispersion, which is quite obvious. At the 
center (r p = 0), the proper motion velocity dispersions have the 
same value, namely 37r(12 — q)/(2‘ (6 — q)). The greater g, the 
higher the dispersions (first derivative is always positive). This 
feature was also noticed in paper I for the line of sight velocity 
dispersion and for the square of the transverse velocity. 

With these values of the velocity dispersions, the total ki¬ 
netic energy T can easily be calculated. Let <r/(r p ) be the 


sum of the three projected velocity dispersions, so af(r p ) = 
cr p (r p ) + o"t(r p ) + cr~(r p ), and cr~(r) be the sum of the three 
spherical velocity dispersions given by 

w( r ) = + r 2 )“= (41) 

and 

4(' r ) = ^(r) = (! + r 2 )' 7 (1 - f ). (42) 

so <r 2 s (r) = <Tr(r) + <r||r) + a |(r), then: 

y + oo 

T = 7r / P p (r p )<rf(r p )r p dr p 

Jo 

= 2tt j p(r) a 2 b (r) r 2 dr = . (43) 

The kinetic energy is independent of g, by virtue of the virial 
theorem. 

Radial cluster types have their greatest contribution to the 
kinetic energy at the center, considering the previous remark. 
Since the total kinetic energy is independent of g, tangential 
cluster types must have their prominent contribution at the 
edge. This feature is affirmed in the plots of the velocity dis¬ 
persions in Fig. 5. 

Moreover, a simple linear relation exists between the three 
velocity dispersions, namely 

o-l(rp) = o- 2 t (r p ) + 12 q _ q O-l(r P ). (44) 

This relation is independent of the projected distance r p , and 
gives again a distinction between both cluster types. For radial 
clusters the line of sight velocity dispersion is always greater 
than the transverse velocity dispersion varying from equality 
if q = 0 until cr p = a"j: + <r”/5 if q = 2. On the contrary, in 
tangential clusters the transverse velocity dispersion is always 







Fig. 6. The dependence of c /( and 7 1 as a function of r p for various q values. 


greater than the radial velocity dispersion, now varying from 
equality in the isotropic case until a 2 = a p + a 2 for q — 00 . 

When we compare the projected velocity dispersions with 
the spatial velocity dispersions, given in Eq.(41) and Eq. (42), 
we notice the similarity between cr~ ( r p ) and a ” (r) on the one 
hand, and the dispersions cr 2 (r p ) or cr 2 (r p ) and cr p (r p ) or <r|(r p ) 
on the other hand. So we define the parameters 7 1 and -) p 
analogous to Jf (cfr. Eq. ( 6 )): 


7t = 1 - 


a t(r P ) 
<j 2 c ( r p ) 


7p = 1 - 



g 5 4 = , 1 . 12 ' 5 4 

12-q l + r p 12 —q 1 + rl' 

q 4 r P ~ 1 = (_i + 12 ) 4 r P ~ 1 

12-q 1 + r 2 12 -q’ 1 + r 2 


(45) 
• (46) 


In Fig .6 these relations are plotted for various q values. The 
parameter 7 1 is a good diagnostic for the orbital structure: its 
sign is the same as the sign of q. In the center, however, it is 
zero for every value of q. The parameter -) p on the other hand 
is in the center region (r p < 0.5) positive for negative q values 
and negative for positive q values, and is zero for r p = 0.5. 
Therefore it is the parameter to use when only dispersions in 
the center are available. For r p —*■ 00 the parameter 7 1 varies 
from —5 for q — 00 to 1 for q = 2 and the parameter -) p 
varies from —4 for q —*■ —00 to 0.8 for q = 2. So at the edge, 7 1 
delivers us a better diagnostic since its range is wider. As usual, 
central velocity dispersions are less useful for determining the 
orbital structure than velocity dispersions obtained at the edge 
(in the unlikely case that both have the same accuracy). This is 
clearly visible in Fig. 6 , where the range in -) p and 7 1 is smaller 
in the center than at large radii. 

From the relation (44) another peculiar result can be de¬ 
duced, namely 


Ip = It — 


g 

12-q 


(47) 


which means that the left picture of Fig .6 can be obtained from 
the right one by pushing the curves up for negative q and down 
for positive q. 


5. The curtosis 

In the same way as for the velocity dispersions, the fourth order 
moments 77 ) (r p ), Tf(r p ) of the marginal distributions can be 
derived. They are: 


168 - 22q + q 2 

~ 2 - 5 - 7 - (7 — g)( 6 - g)(l + r 2 )’ 


(48) 


168 - 22q + q 2 - 4g(34-g)i 2 + l&q(q + 2)x i 
2 - 5 - 7-(7 — g) (6 — g)(l + r 2 ) 


(49) 


whereby x 2 = r p /( 1 + r p ). Instead of the fourth order mo¬ 
ments, one preferably considers the dimensionless curtoses 
cur c = 7"c ( r P )/ cr tA r P ) and curt = r 4 (rpj/cr 4 (r p ), which are a 
measure of the degree to which the marginal distributions are 
peaked: the greater the curtosis, the more peaked the distri¬ 
bution is. Usually this value is compared with the Gaussian 
distribution for which the curtosis is 3. 

A remarkable feature is that the curtosis of the marginal 
distribution in the center direction is independent of the pro¬ 
jected distance r p . Figure 7 shows us the dependence of this 
curtosis on q. The curtosis is in every point smaller than three, 
so less than the normal distribution. Also radial clusters are a 
bit more peaked than tangential ones. 

Figure 8 shows us the curtosis for the marginal distribution 
in the transversal direction. At every projected distance r p , 
radial clusters have a greater curtosis than the tangential ones, 
as can also be noticed in the lower panels of Fig. 4. This was not 
the case for the curtosis of the line profile, which gave us in the 
center region a totally different picture, namely greater values 
are obtained for the tangential clusters! Just like the case of 
the velocity dispersions the scope at higher projected distances 
is greater for the marginal transversal distribution than for the 
line profile. This means again that from the proper motion one 
can get more accurate information. 

Another peculiar feature are the very high values of the 
curtosis in the transversal direction for very radial systems. In 
the limit for r p —*■ 00 or for 1 ^ 1 the transversal curtosis is 
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distance from the center and on the orbital structure. Combin¬ 
ing these three velocity dispersions, it is in principle possible 
to determine the distance. 

For the Plummer sphere, the distance can be obtained very 
easily as follows. Let fi c and fit be the proper motion compo¬ 
nents as defined in sec. 2.2. Since 7 t and 7 P are defined as a 
function of projected velocity dispersions, they are observable 
at each distance r p . Moreover 7 1 is independent of the distance 
to the system A, so the relation (45) provides us with a q for 
each r p . The values of q obtained this way, will be independent 
of r p only in the Plummer case. With this value of q , -) p can 
be calculated, and since 


q 

Fig. 7. The curt.osis of the marginal distribution in the center di¬ 
rection as a function of q. 



Fig. 8. The curt.osis of the marginal distribution in the transversal 
direction as a function of the projected distance. Different values for 
q are indicated. 


2 11 (6 — g) (4 — q)/( 7 r 2 • 3 • 5 • 9(2 — q)(7—q)), which means that for 
q —s 2 this expression tends to infinity. So the two dimensional 
distribution becomes more and more one dimensional along 
the center axis, and has a very high top value. This has been 
checked by calculating the distribution (by numerical integra¬ 
tion) for r p = 100 and q = 1.9. The shape of the distribution 
for large r p and q close to 2 is another confirmation of the ear¬ 
lier mentioned bad behaviour of the series approximation for 
these q values,. 

6. The determination of the distance to a Plummer 
sphere 

The proper motion dispersions (to the center of the system 
and perpendicular to it) depend on the distance to the system, 
on the projected distance from the center and on the orbital 
structure (q). The line of sight velocity dispersion however is 


> = 1 - dr = 1 - 


A 2 fic 


(50) 


with (7 p known and independent of A, the distance A can be 
determined. 

The relative error 8 A/A can be found by differentiating 
the first, of the Eqs. (45) and (46) towards the dispersions and 
the second towards the parameter q. The result is 


— _ —T _|_ rjp 1 - 7t I ‘‘’fit + | Ip — It | 8tic 
A <j p 7t 1 — 7p fit 7t( 1 — Ip) fic 

or expressed in measured quantities: 


(51) 


8A _ 8a p | | 4fp — 1 A 2 //,t | &fi t 

A <7 p 5 r p a p fi t 

+ l-l + ^f^^l —, 

5 r p (7 p fi c 

which is formally independent of q. 

Equation (51) will not hold when 7 t = 0. This occurs either 
when r p = 0 or q = 0. In the former case (the center of the 
system), no appropriate measurements can be made* Even Eq. 
(52) shows this effect. The function of r p in this expression 
varies from 0.6 for r p = 1 until 0.8 for r p going to infinity. 
The greatest influence on the relative error is then due to the 
factor A 2 /(t/(Tp, which makes it very clear that the larger the 
distance A, the greater the error on the distance. In the latter 
case (isotropic), the three dispersions are all equal to: 


(52) 



(53) 


The distance A can then be obtained from A 2 //, 2 = a p or 
A 2 //, 2 = cr p and has the relative error: 


8 A 

8<jp ^ 8fi c 

or 

8 A 

8<jp ^ 8fit 

A 

<Jp He 


A 

&p Ht 


(54) 


This formula comes in handy for a quick order-of- 
magnitude calculation of the relative error on the distance 
determination if the observational errors are known and the 
modelling errors are not too large. 


7. Conclusions 

In this paper we have explored the relation between the internal 
orbital structure of a spherical star cluster and the observed 
proper motion distributions. This we did for a particular family 




carry over, at least qualitatively, to most real clusters. 

We find that simple inspection of the contour lines of these 
two-dimensional distributions suffices to get valuable informa¬ 
tion on the orbital structure, as is obvious from Figs. 2 and 3: 
the shape of the contour lines gives a direct and clear picture 
of the orbital structure. Moreover, we can get this information 
without the knowledge of the distance. 

We also introduce two observable parameters, jt and ~/ p . 
One of them, jt, is the observational equivalent of Binney’s 
anisotropy parameter /?, and measures directly the amount of 
anisotropy. This is the quantitative equivalent of the statement 
in the previous paragraph. 

As a by-product, we show that the moments of the proper 
motion distributions can be used to calculate these distribu¬ 
tions. This means that it may suffice to calculate only the mo¬ 
ments, and be reasonably confident that the proper motion 
distributions are meaningfully constrained. This will be neces¬ 
sary for modelling real clusters, where we can make no direct 
use of the specific formulae given here. 

Last but not least, we show how to determine our distance 
to the system, if all three projected dispersions are known. This 
can be done, because the proper motion components are dis¬ 
tance dependent and the line profile is distance independent, 
whereby they are all three interconnected through the orbital 
structure. The relative error on the distance depends on the 
relative error of the velocity dispersion of the line profile, and 
on the relative error of the two measured proper motion com¬ 
ponents which are multiplied by an appropiate function of the 
distance. 

In a future paper, we will apply the concepts developed 
here to model globular clusters for which both proper motions 
and radial velocities are available. 

Acknowledgements. We would like to thank the anonymous ref¬ 
eree for a careful reading of the manuscript. 


can be written as (see Dejonghe, 1986): 


CM{F} = “ 2 2 M{g(r)}. 

(2tt) 2 F(1 — i/J) 

In the case of the Plummer model this transform is 

3Ff6 — of aiA-T-B+o £ 

CM{F) = \ q l ——- g -— 2 2 M{g(r)}, 


4tt(2 


tt )2 r(i-f) 


(A3) 


(A4) 


and together with the Mellin transform we have 


f+ CO 


M{g(r)} = 


r(f/?)r(i( g -/?)) 


r 13 1 (1 + r 2 ) 2 dr = 


2T(|g) 


,(A5) 


Equation (Al) then becomes: 


r + co 


F(E,L) 
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3r(6-g)V2 1 t*" T(1=Z)T(§)T(3=Z) _ p 

4(27r) fr(Ig) 2*,- J__ P F(l-f) 


1 

2wi 
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e cl 


da. 


— ZOO 


(A6) 


Since the last integral is nothing else than an inverse Laplace 
transform C~ 1 (a q ~ 5 ) = ^> 4-9 /r(5 — q) and the first one is 
Barnes’ integral (Slater, 1966, p22) Eq. (A6) results in: 


r + co 


F{E,L)y/2iP P d( 


3(5-g) ,,,4-q r (q 1 


87T 2 


9 F , - - • 1 • -r 

21 1 2 ’ 2 ’ ’ p 


(A7) 


8. Appendices 

A. The value at the top 

In this case, the value of the energy and angu lar m omentum 
reduces to E = ip p (ip/ip p — ( 2 ) and L = \j2ip p r p (. The 
Laplace-Mellin transform of the distribution function F(E, L ) 
does not depend on the velocity-component (),and is: 




F(E, L) \/2ip p d( 


f+ca 


CO 


'2 i)p d( 


2wi 


da 


> OL — l CO 
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/• /3+ i co 


x - : I dj3 CM{F}e aE L- 13 . 

d/3 — ioo 

So the integration towards £ can be carried out: 

/•+CO _ 


(Al) 


'2il) p dC, e 
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(\/2^r p () 13 
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(A2) 


To calculate the PM-distribution at the top we still need to 
perform an integration through the cluster, which is 


/• + CO 


^-“dz = 


/, 4 -<1 


r r V( r2 ~ J 

r / 3 —q r(f-f) 

r**„ w ^ ) 


: dr 
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so that Eq. (A7) leads to 


P m r (0,0) = 


3(5 —g) r(f-s 


40F r( 2 -f 


-^2^1 • (A9) 


2 2 


With the relation (Exton, 1978, p!7) 


2 F\(a, b; c; z) = (1 — z) a 2 Ti(a, c — b; c;-) 

1 — z 


(A10) 


the above expression can be transformed into 


pm (0,0) = 


3(5—g) r( 


40F r( 2 -f 


3 _SL\ 

-^(i, 1; — 

1 v o ’ o ’ ’ 


,q i 


(All) 


If the density p(ip, r), which is a function of the potential ip which gives us the value at the top as a simple hypergeometric 

and the distance r, can be expressed as a product of the form function, with argument r 2 /( 1 + r 2 ). 


The moments of the distribution function F(E, L ) are defined 


form the (s + i)th derivative into a contour integral around 
r 2 , so that the summation according to i can be calculated, 
yielding: 


Pn,m,l = M /// F(E, L) v™ dv r dv v dvp 


with M the mass, E the energy and L the absolute value of 
the angular momentum. The integration is carried out over all 
the possible values of the velocities. All odd moments are zero 
because the total system has no net streaming. The anisotropic 
moments can be defined as: 


r(g+|)r(p+|)2 p+s i 
7r r(p + s) 2m 


po,o(ip', t) dip' 


This is a general expression, valuable for all kind of fio ,o (ip, t) ■ 


fl 2 n, 2 j = 2ttM JJ F(E, L ) v 2n t)y +1 dv r dvr (B2) 

with v r the velocity component pointed towards the center 
of the system and vt the tangential velocity : vt = (u/, + 
Vp) 1 ^ 2 . The relation between both defined moments can easily 
be shown by putting Eq. (Bl) in cylindrical coordinates: 


The expression of the mass density in the Plummer model 
is given by fto,o (ip',t) = (3/47r) ip' 5 ~ q ( 1 + t)~ q ^ 2 , so that the 
integration towards ip' in Eq. (B4) is a simple beta-function, 
namely: (/>(r)) p+s+6-q B(% — q, p + s). Let the result of the 
integration through the cluster be represented by (v c vt) 2p,2s , 


1 D . 1 , 1 . 

H2n,2m,2l — ~ B(Ul +—,(+—) H2n,2(m+l) 


The anisotropic moments can be expressed in function of the 
potential and the distance r of the center through (Dejonghe, 
1986): 

H2n,2m{lp, r) 

2 m+n r(ra + |) 
r r(m + ra) 


To calculate this integration, we consider new variables v and 
u by means of v = (t — r/)/( 1 + r/) and u = (r 2 — r/)/( 1 + r/), 
and another quite obvious notation for the product of gamma 
functions (Slater, 1966, p41), so that Eq. (B9) leads to: 

/ \2p,2s 3 • 2 P+ “ r f p+ s+ 6 — q\ ( 

{VcVt) = P + S + 6- q (1 + rp) 


x / (/’-/’ 


—1 nm / 2m 


Dp (r 2rn no,o {ip' , r))) dip' 


/•+CO 

/ «“2(1 + 1l)-2 (P+S + 5 - q) dll 

Jo 


with m + n > 1 and with D™ representing the mth derivative 
towards r 2 . 

We wish to give an expression of the moments of the distri¬ 
bution function of the proper motion. So we have to transform 
the velocity components to the coordinate system as defined 
in section 2.2. We have: 


tT(l + u)-i 


= (av r + bvp) 11 vP = a p ' b 


' h ' v p ' v' v v'p (B5) 


with a = r p /r and b = z/r. Taking the mean value of this 
expression, considering that all odd moments have to be zero, 
so the indices p, i and s need to be even, we have : 


2P \ 2(p-i) ,2i 


i l 2(p-i),2i,2s 


Jc W (»-«) P+ ‘ +1 1 7 

with x 2 = Tp/il+Tp), which is nothing else than the cumulative 
mass density. The special case p = s = 0 yields {v c vt} 0, = 
l/(7r (1 + r/) 2 ) which is nothing else than the projected mass 
density at the location r p . 

It is necessary to normalise the moments towards the max¬ 
imum velocity and towards the projected mass density. Its no¬ 
tation will be {xv) 2p ’ 2s and with the aid of Eqs. (9) and (11) 
this normalisation is defined as 

(XV) 2P ’ 2S = ^(l + ^)' (P+i)+2 (M t ) 2p ' 2i (Bll) 

so that (x 7?)°’° = 1. To calculate the integral in Eq. (BIO) the 
expression (v + x 2 ) s is expanded in YH-o (j) x 2t v s ~ l and the 
contour integral is changed back into the derivatives, yielding: 


With the use of Eqs. (B3) and (B4), this becomes: 
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distribution of the proper motion are given by a finite double 
summation. By changing the summation index j to k with 
k = j — i and so m = p + s — i, the last summation is a 
combination of gamma functions and a finite hypergeometric 
series with three parameters in the numerator and two in the 
denominator. The above expression of the moments becomes : 


( X r,) 2p ’ 2s = 3 r /p+|, s + |, 6-g 

' 40F V p+s+6-q 

x f A ^ 2 i (~ 3 ) AAp + s) + i + 2) (f) 8 

“ vv r(z+i)r(|(p+s+5)+z) 

f + + | , -(p + s-i); 

i + 1 , |(p + s + 5) + i; 
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This is in fact a double summation, one with argument x 2 
and one with argument 1. This expression can be transformed 
into a Kampe de Feriet function by considering the following 
two properties: 

— the relation which transforms a a+iFb+i hypergeometric 
function into an integration of a aFb hypergeometric func¬ 
tion is: 

a+iFb+i (c, (a); d, ( b ); z) 

= T(d) 

T(c)T(d-c) 

x / t c_1 (1 — t ) d ~ c ~ 1 aFb ((a); (b);tz) dt (B14) 

Jo 

with (a) = ai,ci2,...,aA the different factors in the nu¬ 
merator and ( b ) = b\ , 62, • • • , 1>B the different factors in the 
denominator. 

— a finite 3 +2 series with argument 1 can be transformed 
into another 3 +2 series with argument 1 after multiplying 
the series with some appropriate gamma functions (Slater, 
1966, Chap.4). 

If we choose out of all the other possible series, that one 
for which the expressions c and d in Eq. (B14) do not contain 
the summation index i, then it is possible to change the order 
of integration and summation in the expression of {xv) 2p ’ 2s - 
This leads after some calculations to the final expression of the 
moments of the distribution of the proper motion: 
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with F a Kampe de Feriet function (Exton, 1978, p24). This 
is in fact a generalization of the hypergeometric function: it 
is a double summation with argument x 2 concerning the first 
summation index i and argument 1 concerning the second sum¬ 
mation index j. After the symbol F, the numbers of parameters 
in the numerator (upper ones) and in the denominator (lower 
ones) are indicated. The first numbers concern the summa¬ 
tion index i + j, and the corresponding parameters are first 
mentioned. They are separated by a colon from the other two 
groups. The first group concerns the first summation index i 
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a finite summation) and zero parameters in the denominator. 
The second concerns the second summation index j, here with 
two parameters in the numerator, which are 1/2 and — (p + s) 
(so again a finite summation), and one in the denominator, 
which is (p + s — q) /2 + 3. The separation of these parameters 
is now indicated by a semi-colon. Special cases of the Kampe 
de Feriet functions are the four better known Appel functions. 

C. Outer regions of the PM-distribution 

We will put forward an approximation of the PM-distribution 
for velocities close to the escape velocities, this means for the 
outer regions of the PM-distribution. Here is v 2 = x 2 + rj 2 —»■ 1 
so that £ 2 is very small. Since there is a double prescription 
of F(E, L), given in Eqs. (4) and (5), two kinds of approaches 
are needed. 


C.l. Small values of r p 

In the case of small values of r p the greatest contributions to 
the outer regions of the PM-distributions come from the stars 
with a radial velocity close to the escape velocity. For those 
stars the angular momentum vector is small and consequently 
L 2 « E, and we can restrict the hypergeometric function in 
Eq.(4) to the first term: 
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By using the usual symbols and normalisations, the PM- 
distribution becomes: 
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The integration towards £ is a beta-function, so that the above 
expression becomes: 
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Transforming the r- variable into ip considering that 1 fr(, al = 
4 /i max 4/(/+ ^ )> leads to. 
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Integrating from 0 to 1 with the substitution ip = ip p (v 2 -\~u( 1 — 
v 2 )) results in: 

P m rAx,v) 
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so that the first factor in the denominator goes to 1 and the last 
one goes to \pl. The integration towards u is a beta-function, 
namely B(j, 5 — q) so that we finally have: 
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C.2. Large values of r p 


In this case the most important contributions to the outer 
regions of the PM-distribution come from those stars whose 
angular momentum is great. So here is L 2 >> E and so we 
restrict the distribution function in Eq. (5) to the first term of 
the hypergeometric function: 
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since in these regions the 2 -values are very small with regard 
to r p . Also v 2 = x + g is nearly 1. This means that Cf 2 is 
very small. If g is small, then the angular momentum is small 
and we are back in the previous case. If g is large, we can then 
approximate L by: 

L~ q « (2f p y q/2 r~ q (g 2 )- q/2 (C9) 


so that the PM-distribution becomes: 
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Integration towards ^ is a simple beta-function, yielding: 
P m r p (x, v) 
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Transforming the r -variable into ip considering that 1 + r^ nax = 
l/ip p n 4 and integrating from 0 to 1 with the substitution ip = 
ipp(v 2 + u(l — v 2 )) results in: 

P m r p (X, V) 

_ 3r(6-g) r p q (g 2 )-? ip~ q 
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For v 2 —> 1 this leads to: 
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P m r p (X, V) 
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